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ABSTRACT 



Context. Main Belt comets represent a recently discovered class of objects. They are quite intriguing because, while having a Tisserand 
invariant value higher than 3, are showing cometary activity. 

Aims. We study the activity of the Main Belt comets making the assumption that they are icy-bodies and that the activity has been 
triggered by an impact. We try to determine the characteristics of this activity. We also try to determine if the nowadays impact rate in 
the Main Asteroid Belt is compatible with the hypothesis of an activity triggered by a recent impact. 

Methods. Due to the fact that the Main Belt comets can be considered as a kind of comets, we apply a thermal evolution model 
developed for icy bodies in order to simulate their activity. We also apply a model to derive the impact rate, with respect to the size of 
the impactor, in the Main Belt. 

Results. We demonstrate that a stable activity can result from a recent impact, able to expose ice-rich layers, and that the impact rate 
in the Main Belt is compatible with this explanation. 

Key words. Minor planets, asteroids: general - Comets, general - Comets, individual: P/2005 Ul (Read) 



1. Main Belt comets: a new class of objects 

The Main Belt Comets (MBCs) are a recently discovered class 
of objects, composed until now by only four objects (Table [TJ. 
These objects are puzzling because they are orbiting in the Main 
Belt on stable orbits, with Tisserand invariant with respect to 
Jupiter higher than 3, but at the same time they are showing a 
more or less stable cometary activity. Four known objects im- 
plies many more presently inactive or faintly active and as a 
consequence undetectable. It is probable that as a result of the 
surveys being conducted (Gilbert & Wigert 2009) more objects 
will be added to this new class in a near future. 
The first of the MBCs that has been discovered is 133P/Elst- 
Pizarro. At the moment of the discovery, in 1992 (Elst et al. 
1996 ), it has been designated as 7968 Elst-Pizarro. When in July 
1996 the body was found showing a dust tail, the activity was 
initially attributed to an impact (Boehnhardt et al. 1998, Toth 
2000), but when later on the object was found periodically ac- 
tive at the successive perihelions (2002 and then 2007) it became 
clear that a different explanation had to be found (Hsieh et al. 
120041 . P/2005 U l (Rea d) and 176P/ LINEAR followed in a short 
time (Read et al. 120051 Green 2006 ), both active and orbiting on 
similar orbits (Table [TJ- More recently, P/2008 Rl (Garradd) has 
been discovered (Garradd et al. 120081 ). 

The Tisserand invariant values rule out the possibility that these 
bodies are newcomers in this part of the Solar System. The or- 
bits of two of them, 133P/Elst-Pizarro and 176P/LINEAR, can 
be linked with the 2 Gyr-old Themis dynamical family (Zappala 
et al. 1 19901 ) and the smaller and younger (10 Myr) Beagle family 
(Nesvorny et al. [2008b . The orbit of P/2005 Ul (Read) is falling 
just outside the Themis family, and P/2008 Rl (Garradd) appears 



to be a dynamically unstable object, probably formed in a differ- 
ent part of the Main Belt (Jewitt et al. 120091) . 
As for the causes of dust ejection, it became soon evident that 
it is due to sublimation. Not only this activity is recurrent, as in 
the case of 133P/Elst-Pizarro, but it is also lasting for weeks and 
months. Moreover, dust velocities, estimated from the sunward 
extent of the coma, are greatly in excess of what would be ex- 
pected from electrostatic levitation or rotational ejection (Hsieh 
& Jewitt p2006l 

If these bodies formed in-situ, the question of how can they 
still contain ice arises. In the protoplanetary disk, opacity due 
to grains inhibited radiative cooling and raised the mid-plane ki- 
netic temperature (Jewitt et al. 2007). Growth and migration of 
solids in the disk would have caused the opacity to change, so 
moving the snow line, that for a fraction of a million years may 
have been within the Main Belt. In this way, asteroids orbiting in 
the Main Belt could have incorporated water ice upon formation. 
It is well know that many asteroids lack the hydration feature in 
their spectra. This has been interpreted as a sign that their ice 
never reached the liquid phase and so could still be found frozen 
in the interior (Jones et al. Q990; Carvano et al. 120031) . These 
bodies would be possible candidates for the new class of MBCs. 
The recent detection of water ice on the surface of 24 Themis 
(Campins et al. 120101 Rivkin & Emery 120101 ) came as no sur- 
prise, and give a strong support to the sublimation-driven expla- 
nation for the MBCs activity. At the end of 2009, two indepen- 
dent teams observing 24 Themis in the infrared range detected 
the presence of water ice widespread on the surface of the body. 
It has been interpreted as frost, replenished by some subsurface 
source. This detection, and those of Main Belt comets, are con- 
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tributing to definitively shatter the old, traditional view of the 
Main Belt as the realm of rocky bodies as opposed to the icy, 
active bodies that can be found farther. 

The implications deriving from the existence of a third source 
(besides the Oort Cloud and the Edgeworth-Kuiper Belt) of 
"cometary" bodies are remarkable. The water contained in the 
terrestrial oceans is thought to come from impacting icy bodies, 
but the few D/H ratios measured until now in comets have double 
values with respect to the average ratio obtained from terrestrial 
water. A third potential source for Earth water, different from 
long-period and short-period comets, could solve the problem, 
if the D/H ratio in MBCs will prove to be more similar to the 
terrestrial one. 

2. The activity of Main Belt Comets 

Dust ejection due to sublimation implies the existence of both 
ice and dust particles. This leads us to imagine that MBCs are 
comet-like bodies, that means intimate mixtures of ice and dust, 
and not rubble-pile rocky objects with some ice filling the inter- 
stices. If this assumption is coupled with the further assumption 
that the activity is driven by sublimation, then we can simulate 
and study the activity of MBCs through the thermal modeling 
codes used for comets, KBOs and icy bodies in general. 
Prialnik and Rosenberg (2009 ) applied their comet nuclei ther- 
mal evolution model to 177P/Elst-Pizarro. In their paper the au- 
thors demonstrate that deep-buried ice could have survived since 
the Main Belt formation time, and that this ice is most probably 
composed, at least nowadays, of water crystalline ice. It would 
have been almost impossible, for ices sublimating at tempera- 
tures well under 130 K, to survive even under an insulating man- 
tle. Schorghofer too (2008) in his paper assumes that the ice, in 
order to survive for such a long time, must have been buried un- 
der an insulating mantle, and gives estimations of the thickness 
of this dust layer. The same author argues against the possibility 
that the ice be mixed with rocky material rather than dusty ma- 
terial. Rocky surfaces are seldom able to retain ice, due to the 
larger thermal conductivity and the larger molecular free path of 
rocks with respect to dust grains. 

Being obvious that exposed ice cannot last for a long time at 
this distance from the Sun , we make the hypothesis that buried 
ice has been exposed in recent times. The most probable expla- 
nation is that of an impact excavating a crater in a dust mantle. 
This hypothesis, and if the present-days impact frequency in the 
Main Belt is compatible with it, will be discussed in a following 
section. 

Starting from the assumptions stated above, we apply our ther- 
mal evolution model to simulate MBCs activity in the hypothesis 
that an impact has recently happened and that a crater has been 
excavated in the dust mantle. The effect of the impact could have 
been to directly expose a fresh ice-rich layer or, if the crater is 
not deep enough, bring an ice-rich layer closer to the surface, 
making it reachable by the the diurnal/seasonal heat wave. The 
aim of the simulation is to investigate the activity of MBCs in a 
more general way and under different hypotheses, and study how 
long can the activity last, on which parameters it depends, and 
which are the characteristics of this activity when a thin mantle 
is still covering an active layer. 

2. 1 . The thermal model 

The thermal evolution model that has been applied to simulate 
the activity of MBCs is described in this subsection (Capria et al. 
I2000at Capria et al. l2000bl Capria et al. l200Tl De Sanctis et al. 



2003, 2005 ). In the next subsection, the simulation of dust com- 
ponent, dust flux and mantle structure are explained in greater 
detail, because these features have been recently improved in 
this continuously evolving code. 

The model is one-dimensional. The spherical nucleus is porous 
and composed of ices (water and up to two different species, typ- 
ically CO and CCb) and a refractory component. Water ice can 
be initially amorphous, and in this case a fraction of gases such 
as CO can be trapped in the amorphous matrix and released dur- 
ing the transition to crystalline phase. In this work, anyway, we 
will deal only with water ice in its crystalline status. Together 
with Prialnik and Rosenberg ( I2009I ), we deem extremely im- 
probable, for ices sublimating at temperatures well under 130 
K, to have survived till now, even under an insulating mantle. 
Energy and mass conservation is expressed by the following sys- 
tem of coupled differential equations, solved for the whole nu- 
cleus: 



pc^ = V[K ■ VT] + Qt 



(1) 



RT dt 



V[G,- • VP,] + Qn 



where p, c and K are the bulk density, specific heat and thermal 
conductivity, is the energy exchanged by the solid matrix in 
the sublimation and recondensation of the water ice, R is the gas 
constant, P, the gas partial pressure, G, its diffusion coefficient, 
and Q/j is the gas source term due to sublimation-recondensation 
processes. 

Gas diffusion coefficients are computed on the basis of the mean 
free path of the molecules in the pore system; the model ac- 
counts for three different diffusion regimes: Knudsen, viscous 
and a transition one. 

The temperature on the surface is obtained by a balance between 
the solar input and the energy re-emitted in the infrared, con- 
ducted in the interior and used to sublimate ices present on the 
surface. The following expression is used: 



S(l -A s )cos6» dT 

d2 = co-T] + K(T S ) — 

Rr n dr 

h ' -'" 

fLo ' H(Ts)h 2 o ■ e(T s )h 2 o 



(2) 



where S is the solar constant, /?/, the heliocentric distance of the 
comet, A j the Bond bolometric albedo of the surface, e the in- 
frared emissivity of surface, <x the Stefan-Boltzmann constant, 
K the thermal conductivity, and 6 the cometary latitude. 
In the last term of right member of Eq. (O, f^ represents the 
surface fraction covered by water ice, H is the latent heat of sub- 
limation and s is the sublimation rate. 

The surface boundary condition used to solve diffusion equation 
when ice is present on the surface is: 



p = p m '(T) 



(3) 



When no ice is present, surface boundary condition is obtained 
assuming free sublimation at comet nucleus surface: 



P = 



(4) 



When the temperature rises, ices start to sublimate, beginning 
from the more volatile ones; even if the nucleus is initially 
homogeneous a differentiated, layered structure is always 
obtained, in which the boundary between different layers is a 
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Name 


a(ALT£ 


e b 


i(% 


P (year£ 


TP 


133P/Elst-Pizarro 


3.156 


0.165 


1.39 


5.62 


3.184 


P/2005 Ul (READ) 


3.196 


0.192 


0.24 


5.63 


3.166 


176P/LINEAR 


3.165 


0.253 


1.27 


5.71 


3.153 


P/2008 Rl (Garradd) 


2.726 


0.342 


15.9 


4.43 


3.216 



" Semimajor axis 

b Eccentricity 

Inclination 

d Orbital period 

' Tisserand parameter 



Table 2. Properties of the dust distributions 



Property 


Silicatic 


CHON/organic ] 


Size classes 


6 


6 


Size ranges (m) 


10~ 6 - 10-' 


io- 6 - io-» , 


Density (kg m ~ 3 ) 


2.5 • 10 3 


1.0 • 10 3 | 


Conductivity (grain) (K kg -1 nr 3 ) 


3 


0.25 


Conductivity (mantle) (K kg~' irr 3 ) 


0.01 


0.2 



sublimation front. Surface erosion due to ice sublimation and 
particles ejection is taken into account. 

The dust component is represented as size distributions of 
spherical grains characterized by a density, a specific heat 
and a thermal conductivity; in the code we can consider up to 
two distributions, with a maximum of ten size classes each. 
Dust distributions characteristics and their physical properties 
were chosen to simulate as much as possible Giotto findings. 
These experiments onboard of the mission demonstrated that in 
P/Halley's refractory component it is possible to evidentiate two 
chemical phases: an organic one with low atomic weight, called 
CHON, and a Mg-rich silicate one (Jessberger and Kissel 1991 ; 
McDonnell et al. |199U Mumma[T997). Measured particles were 
in the range of 10~ 16 - 10~ n g. One third of it did not contain 
organic components, and the remainder was a mixture of the two 
phases; CHON component was probably coating silicate cores. 
Average densities were found to be 2500 kg rrT 3 for silicate 
grains and 1000 kg irr 3 for CHON-dominated particles. 
For these reasons, the dust component is simulated by two dis- 
tributions: the first distribution is composed of silicatic particles, 
and the second one by mixed silicatic/CHON particles. The 
characteristics attributed to the two distributions are collected in 
table [2] and are explained hereafter. 

These two phases have different physical properties (density, 
thermal conductivity, and so on); as a consequence, the charac- 
teristics (thermal conductivity, tensile strength and cohesion) of 
a crust mainly composed of silicatic particles would be different 
from those of a crust with an abundant component of refractory 
organic material. 

For the grains of the first distribution a thermal conduc- 
tivity value of 3 W Kr'rrT 1 , typical of silicates, has been 
adopted, while for the CHON-dominated grains a value of 0.25 
W Kr'irT 1 , typical of organic substances thought to be ana- 
logues of cometary organic material (Komle et al. 1996), has 
been preferred. Average density is 2500 kg irT 3 in the first dis- 
tribution, and 1000 kg m -3 in the second one, accordingly with 
Giotto findings. 

The size distributions of dust particles in the nucleus could be 
different from those measured in the coma: observed particles 



are the smaller ones, and larger particles can be seen only by 
radar (Harmon et al. I1999I I. The observed exponential distribu- 
tions could be the results of fragmentation processes in the inner 
coma. A way to define an initial size distribution in the nucleus 
could be to suppose that the grains are the result of the coagu- 
lation processes acting in the protosolar nebula: following this 
idea, we attributed to both kind of grains a gaussian distribution 
taken from a paper by Coradini et al. ( 1 19771 1 on the accretion of 
grains in the solar nebula. 

It remains to define a value for the conductivity of a silicate- 
dominated crust and for the conductivity of a crust with a high 
abundance of organic refractory material. It was found, during 
KOSI experiments (Grim et al. 119931 1, that thermal conductivity 
of dust layers composed of silicatic grains is very low (~ 10~ 2 - 
10~ 3 W Kr'rrT 1 ): even a thin layer of dust causes a large tem- 
perature gradient. Following this result, we attributed a value of 
0.01 W K -1 irr 1 to the thermal conductivity of a crust composed 
only of silicatic particles. In the case of the "organic" distribu- 
tion, lacking for direct measurements, we deduced the physical 
properties from the results of the laboratory experiments con- 
ducted at the Space Research Institute of Graz (Komle et al. 
1996). The thermal evolution of an analogue of cometary ma- 
terial composed of a mixture of minerals, ices and hydrocarbons 
has been studied, obtaining a sort of cohesive crust composed 
of minerals glued together by organics. This kind of crust has a 
thermal conductivity higher than that of silicatic crust, and also 
has a greater cohesive strength: an admixture with organic mate- 
rial raises crust conductivity of one order of magnitude or more 
with respect to a crust composed of silicatic grains. So, our "or- 
ganic" crust is characterized by a higher thermal conductivity 
than the "silicatic" one: we attributed to a crust mainly com- 
posed of organic particles a value of 0.2 W K^'irT 1 . 
At the beginning of computations the dust particles are embed- 
ded in the porous ice. When the ice sublimates the embedded 
particles become free, in a number proportional to the amount 
of sublimated ice, and can undergo the drag exerted by the gas 
flux. Average pore size is increasing due to the ice sublimation: a 
"pore widening factor" proportional to the amount of sublimated 
ice is locally applied to the pore average size. In this way we are 
taking into account that a sublimating ice produces holes and 
cavities, and as a consequence free particles can move towards 
the surface (particles can move and be blown off only if their 
size is less than the local average pore size!). Free particles that 
are found close to the surface can move toward the surface and 
be blown off or accumulate to form a dust mantle. Free particles 
are moving under the opposite effects of gravity and gas flux; 
only grains smaller than the pore average size can move toward 
the surface in the pore system. When the ice is condensing, a 
proportional amount of free particles become embedded and the 
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pores becomes smaller (porosity decreases). 
To determine how many particles can be blown off to contribute 
to dust flux, and how many can instead be accumulated on the 
surface, the different forces acting on the single grain are com- 
pared, obtaining for each distribution a critical radius, depending 
on grain density. The critical radius represents the radius of the 
largest particle that can leave the surface of the comet: 

. 3 ®h 2 o ■ Vh 2 o r , 

a ~ A M L»lJ (5) 

4 PdKif[G-^ -R n co 2 cos 2 8] 

where G is the gravitational constant, M„ the mass of the comet, 
R„ its radius, a> its angular rotation velocity, pd us , the dust grain 
density, <1>h 2 o is the gas flux, and Vh 2 o the velocity of the gas 
flux. In the Eq. © the numerator represents the lifting force ex- 
erted by the outflowing gases, and the denominator represents 
the gravitational attraction corrected by the centrifugal force. At 
each time step the number of free dust particles for each size 
class and the value of the critical radius are computed. All the 
free particles with radius a < a* are considered ejected (blown 
off) and contribute to the dust flux, while those with a > a* 
accumulate on the surface. This refractory surface layer (dust 
mantle) is not necessarily stable, and can be later on blown off if 
the gas flux is strong enough. It could also increase in thickness 
if more and more particles are not able to leave the surface, sub- 
stantially reducing and then quenching the gas flux. Dust mantle 
compaction due to a volume reduction when ice is sublimated is 
taken into account. If the dust mantle is not too thick the gas flux 
from the interior is able to entrain a dust flux, so it is possible to 
have a dust flux even when there is no ice left on the surface. 
In this way, even an originally homogeneous body becomes 
quickly highly inhomogeneous: composition, density and poros- 
ity are changing with depth. 

2.2. Thermal modeling of P/Read 

A thermal model depends not only on the assumptions of the 
model itself, but also on the values attributed to the initial 
parameters. These parameters describe the orbit, the initial 
state of the body and the properties of the matter of which it 
is composed. They are defined or derived, when possible, from 
observation and laboratory experiments. A given set of these 
parameters define a "Case". By changing a small number of 
these parameters, and keeping fixed all the rest, we can build 
different Cases that are the subject of the simulation. 
To apply the above described model to a real case we need to 
define the characteristics of the body that we want to simulate, 
assuming an initial structure and composition. As for the real 
case, we chose to simulate the behaviour of P/2005 Ul (Read). 
The results of this simulation can be applied also to the other 
MBCs known until now, because orbits and sizes are reasonably 
similar. 

As written before, we start from the assumption that the model 
body has a cometary nature, that is it is composed of ice 
and dust grains. A further assumption is being made that the 
surface of the body is covered by a devolatilized mantle, and 
that a recently happened impact has opened a crater in this 
mantles. As a consequence, the heat from the Sun is now able 
to reach ice-rich layers, so triggering a cometary-like activity 
(sublimating gas entraining dust particles). At the beginning of 
the simulation, the immediate effects of the impact (release of 
heat) have already ceased. 

As for the composition and physical properties, it is difficult 



Table 3. I. Initial value of the parameters defining the model 
body. 



a 


3.165 


e 


0.253 


Dust/ice 


3 


Initial temperature (K) 


130 


Surface albedo 


0.04 


Surface emissivity 


0.96 


Spin period (hours) 


10 


Initial pore radius (m) 


10" 4 


Initial porosity 


0.4 


Initial radius (km) 


0.3 



to make safe assumptions. We have no reason to think, any- 
way, that these bodies are absolutely similar to the classical 
comets (Jupiter-family, Halley-family and Long Period comets), 
because their origin and evolution history is different. They 
formed closer to the Sun than the comets now stocked in the 
Kuiper Belt or Oort cloud. During the Solar System formation 
period materials from different distances from the Sun were 
mixing up and the snow line could have been closer to the 
Sun than nowadays: this can explain the presence of dust-ice 
bodies orbiting what is now the Main Belt. Collision evolution 
also played a major role in the past of these bodies, probable 
survivors of the catastrophic event that gave origin to Themis 
family. Hardly the ice now remaining in MBCs was formed 
in the amorphous state, because this requires extremely low 
temperatures (much lower than 50 K). For similar reasons, it 
is extremely improbable that this ice could still be mixed with 
other ices more volatile than water that, even if they existed, 
should have been disappeared long time ago. We will make the 
assumption that the model body contains only H2O crystalline 
ice. 

Some of the values given to the initial parameters of the model 
are listed in the table [3] The properties attributed to the dust 
distribution are listed in the first column of the table [2] The 
values given to orbital parameters (semimajor axis and eccen- 
tricity) and to the radius of the model body are derived from 
observation, and refers to P/2005 Ul (Read). The values given 
to surface albedo, emissivity, and pore radius can be considered 
typical of comets, and are largely used in comet modeling 
(Huebner et al. 2006). The spin period is not known, so we are 
using a value (10 hours) typical of many small bodies. The value 
for dust/ice ratio we are using in this simulation is 3, higher than 
the one usually assumed for comets (1), as a way to take into 
account the origin of the comet in a zone of the Solar System 
less rich in volatiles. The porosity is 0.4. The initial temperature, 
130 K, has been defined taking into account the orbit of the body. 

The dust component is described as a distribution of 6 size 
classes with similar density and thermal conductivity (see table 
0. 

The model was run in the fixed-latitude approximation. The 
reasons for this choice are that the rotation period is not known, 
and the variations of gas and dust production with respect to the 
daily hour are irrelevant for the present analysis. In this way, gas 
and dust production rates depend only on the latitude and on the 
position along the orbit. 
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Fig. 1. Case 1 - Water flux along the orbit for a P/Read-like body 
on which an impact has exposed ice-rich layers. Upper left panel: 
impact at the equator, water flux with respect to time, for the first 
140 years after the impact; upper right panel: impact at the equa- 
tor, water flux with respect to time for a period of one hundred 
years immediately before the disintegration of the available icy 
material; bottom left panel: impact at the equator (upper curve) 
and at 30° (lower curve), water flux with respect to heliocentric 
distance. The fluxes are shown 10 years after the impact; bottom 
right panel: impact at 85°, water flux with respect to heliocentric 
distance. The flux is shown 10 years after the impact. 



2.3. Results: exposed ice 

Case 1 has been ran assuming that the impact has directly ex- 
posed to the sun heating an icy zone located at given cometo- 
centric latitude. The model is ran for 4 different latitudes of the 
impact crater: equator (sub-solar point), 30°, 60° and 85°. The 
dust component is represented by one distribution, the silicatic 
one (first column of table |2j. With these parameters, bulk ther- 
mal conductivity along the whole evolution has an average value 
of 0.07 K kg -1 irT 3 ; average density values for the matter 
composing the model body are of the order of 586 kg m -3 . 
The results can be seen in the figures [1] [2] [3] and |4] showing, 
respectively, gas and dust fluxes, stratigraphy and internal tem- 
perature profile. After a time of almost 900 years, the ice ma- 
terial exposed by the impact at the subsolar point is completely 
consumed, and gas flux, that has been continuously reducing (Q~| 
upper panels), is completely quenched. At all the latitudes, most 
of the dusty material is left behind: not all the dust grains freed 
from the ice are going to contribute to the dust flux. At higher 
latitudes the situation is different, as can be seen in the bottom 
left panel of figure [1] showing the gas flux with respect to helio- 
centric distance for impacts at the equator, 30 deg and at 60 deg. 
In the bottom right panel of the same figure the water gas flux 
with respect to heliocentric distance is shown for an impact at 
a latitude of 85 deg. In this case, the flux is extremely reduced. 
Bottom panels show the gas flux ten years after the impact ex- 
posing the ice-rich layers. 

The dust flux with respect to time for an impact at the equator on 
a period of 60 years, starting 10 years after the impact, is shown 
in figure [2] The reduction of the dust flux from orbit to orbit can 
be clearly seen. For impacts at higher latitudes, the dust flux is 
correspondingly reduced, similarly to what happens to the gas 
flux, and is almost null when the impact crater is supposed to be 
located close to the pole. 



Fig. 2. Case 1 - Impact at the equator, dust flux with respect to 
time for a period of 60 years following the impact. 

The erosion level of the surface of the crater varies with 
the latitude of the impact and is proportional to the intensity 
of dust and gas emissions, as can be seen in figure [3] showing 
the stratigraphy (surface radius, in meters, with respect to time 
in years) for a period of 50 years after the impact: it is almost 
null close to the poles (lower panel), where the solar input 
is minimum, while at the sub-solar point the available ice is 
consumed in less than one thousand years (almost two meters of 
surface lost per orbit, upper panel). 

The temperature profiles shown in figure [4] are photographing 
the situation at the aphelion, after almost 9 hundreds years of 
revolutions around the Sun. At the equator only 30 meters of 
active material remain (panel a), while at 85 deg of latitude 
(panel b) not even 1 centimeter of surface has been lost. Surface 
and interior temperatures give an explanation for this difference: 
at temperatures equal or less than 145 K, as is the case at 
latitudes close to the pole, water ice can be considered stable 
against sublimation. At equator, in the 30 meters of icy material 
remaining after 900 years of activity, interior temperatures are 
much higher, always well above the sublimation limit. The first 
5 meters under the surface are under the influence of the orbital 
heat wave. In the deeper layers the temperature slowly decreases 
until it reaches 155 K, 25 K higher than the temperature at the 
time of the impact (130 K). 



2.4. Results: mantle-covered ice 

In order to trigger some activity, it is not necessary that the 
impact exposes directly icy-rich layers to the action of the Sun. 
The result of a shallow impact in a devolatilized mantle could 
be excavate a crater, without directly exposing ice-rich layers, 
but leaving only a thinner dust layer covering regions that had 
remained safely buried for millions of years. In this case it could 
be possible not only activate the body (as a consequence of 
the impact itself) but also generate a stable or periodic activity, 
although fainter than in the situation described in the preceding 
section. It is well known that ice sublimation can be possible 
from under a mantled surface (see, for example, Capria et al. 
120011) . providing that ice-rich layers are reasonably close to 
the surface and the mantle itself is somewhat porous. This 
phenomenon has been observed, for example, on the comet 
9P/Tempel 1, target of Deep Impact. The few small patches 
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Fig. 3. Case 1 - Stratigraphy showing the erosion of the sur- 
face with respect to time. Upper panel: impact at the equator 
(lower black curve) and at 30°(red curve); lower panel: impact 
at 60 "(lower black curve) and 85 "(orange curve). At the begin- 
ning of computations, the surface has always a radius of 300 m. 
Depending on the latitude of the impact, the radius shrinks with 
respect to time. 
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Fig. 4. Case 1 - Internal temperature profile for a P/Read-like 
body at the aphelion, 900 years after the impact; impact at the 
equator (panel a) and at 85° (panel b. 



of icy surface visible in the returned images are not enough to 
account for the activity of the comet, if sublimation from the 
mantled zones is not considered (Belton and Melosh 2009). 
To simulate this kind of activity, 4 Cases have been ran with a 
different assumption than in the preceding section: the bottom 
of the crater is still covered by a dust layer that has survived the 
impact. In this case, the duration and intensity of the activity 
depend strongly on the mantle thickness and physical properties. 
It would be impossible to run Cases exploring all the different 
combination of input parameters having a direct influence on 
the results. Two different thicknesses (0.1 and 0.5 m) and two 
extremely different compositions (silicatic and organic/CHON) 
have been considered for the surviving dust layer, obtaining 
4 Cases (see table [4] The thickness of the mantle before the 
impact should have been much larger: in order to keep sealed 
volatile materials on an extremely long timescale, tens of meters 



of insulating materials are necessary (see also Schorghofer 
2008). The surviving dust layer thicknesses to be tested have 
been chosen somewhat arbitrarily, on the basis of our preceding 
experience on comet nucleus modeling. The results of these 
Cases can be easily extrapolated to different thicknesses of the 
surviving dust layer. 

The aim of testing two different dust compositions is instead 
to show how the characteristics of the activity would depend 
on the physical properties of the involved materials: we are 
aware that pure silicatic or pure organic dust are somewhat 
unrealistic. The different compositions correspond in the code 
to two dust distributions with the characteristics listed in table 
[2] The remaining input parameters are the same as the ones 
listed in [3] The impact point is always assumed at the equator 
(sub-solar point). 

Figures [5] |6] and [TJhow the result of these 4 Cases. In fig. [5] 
the gas flux is shown for a body with a 0. 1 m thick surviving 
mantle layer, assuming a silicatic composition (upper panels, 
Case 2) and an organic composition (bottom panels, Case 4). 
The left panels show the evolution of the water flux in the years 
following the impact, while the right panels show the flux after 
almost 500 years after the impact. A similar figure © show the 
result of Cases 3 and 5, where the thickness of the surviving 
mantle is 0.5 meters. 

A quite obvious result is that the water ice is sublimating with 
a rate much lower than in Case 1, depending on the thickness 
of the surviving dust layer: between Cases 2, 4 (0.1 m thick 
layer) and Cases 3, 5 (0.5 m thick layer) there are orders 
of magnitude of difference. Less obvious are the long-term 
evolution, and the dependence of the activity on the properties 
of the material forming the dust layer. When the computations 
begin, the interior temperature of the body is homogeneous 
and equal to the given initial temperature, 130 K in this case. 
With the passing of the years, the interior temperature slowly 
increases, and this corresponds to an increase in the sublimation 
rate. The process goes on till the interior temperature stabilizes, 
along with the sublimation rate (right panels in figures [5] and 
[6] Surface temperatures, mainly determined by solar input, are 
spanning the range from 148 K to 191 K in all the four Cases. 
Due to the action of two competing processes, the thickness 
of the surviving dust layer changes, even by few millimeters, 
in the time after the impact. After 5000 years, in Cases 2 and 
4 the thickness of the mantle has increased, at the bottom, 
by 15 centimeters and almost 12 centimeters respectively. 
This increase is the result of the slow devolatilization of the 
icy-reach layers adjacent to the bottom of the mantle. During 
the same time, a faint dust flux has eroded the surface by 2 and 
3 millimeters in the two Cases. As for the Cases in which the 
initial thickness of the mantle was higher, after 5000 years of 
evolution less than 5 centimeter of devolatilized material has 
been added at the bottom of the crust, while the surface erosion 
is less than 2 millimeters in Case 5 and nothing in Case 3. 
The different values attributed to the thermal conductivity and 
density of the refractory component of the model body (see 
table|2]i explain the differences in the activity between the Cases 
2,3 and 4, 5. Due to the higher thermal conductivity of silicatic 
grains, in Cases 2 and 4 higher internal temperatures are attained 
with respect to Cases 4 and 5. Figure [7] shows the temperature 
profiles of Cases 3 and 5, at the beginning (left panels) and 
towards the end of computation (right panels). The strong 
temperature gradient between the surface and the sublimation 
front (where water ice, mixed with dust, still survives), is visible, 
along with the difference in surface temperature between the 
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Table 4. Input parameters characterizing the surviving dust layer 
at the point of the impact in the 5 Cases 



Case 


1 


2 


3 


4 


5 


Dust distribution 


silicatic 


silicatic 


silicatic 


organic 


organic 


Thickness of the dust layer at impact location (m) 





0.1 


0.5 


0.1 


0.5 




Fig. 5. Cases 2 and 4 - Water fluxes with respect to time, at the 
beginning and at the end of computation. The thickness of the 
surviving dust layer is 0.1 m. Upper panels: Case 2 (silicatic 
dust). Bottom panels: Case 4 (organic dust). 
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Fig. 6. Cases 3 and 5 - Water fluxes with respect to time, at the 
beginning and at the end of computation. The thickness of the 
surviving dust layer is 0.5 m. Upper panels: Case 3 (silicatic 
dust). Bottom panels: Case 5 (organic dust). 



aphelion, at 3.96 AU, and a point on the orbit, 3.16, closer to 
perihelion (2.36 AU). 

At a certain point, the increase in the internal temperature 
ceases because an equilibrium temperature has been reached. 
This process has a different duration depending on the physical 
properties of the materials forming the body. It can be noticed 
that when the computation stops, 5000 years after the impact, the 
model bodies are still active, and it is reasonable to assume that 
they will continue to be active for a very long time. As for the 
thickness of the surviving dust layer, there are two competitive 



Fig. 7. Cases 3 and 5 - Temperature profiles at equator. Left up- 
per panel: Case 3, 1700 years after the impact, aphelion. Right 
upper panel: Case 3, 5000 years after the impact, 3.16 AU. Left 
bottom panel: Case 5, 1700 years after the impact, aphelion. 
Right bottom panel: Case 5, 5000 years after the impact, 3.16 
AU. 



processes at work: dust emission from the surface (and the clos- 
est layers), and progressive devolatilization of layers close to the 
bottom of the mantle. The second process is the winner, and so 
the thickness of the mantle is slowly growing. As a consequence, 
the gas flux is slowly reducing. In order to study how much time 
would the body remain active, much longer, and heavy, compu- 
tations would be required. This will be the subject of a future 
work. 



2.5. Comparison with observation 

The results of the theoretical models should now be compared 
with observations. Only in this way would it be possible to judge 
which ones, between the different Cases presented in the above 
sections, are closer to reality. Unfortunately, what we have from 
observation about P/Read is not enough to make possible such 
a comparison. From the paper by Hsieh et al. (/2009), who ob- 
served this object for almost three months in 2005 and one night 
in 2007, P/Read has been always found active in 2005, while in 
2007 it was too faint to permit an estimation of the activity level. 
It must be noted that P/Read in July 2005 passed at its perihe- 
lion, and the aphelion passage was in May 2008, a few months 
after the observation of 2007. In 2005, coma and tail were clearly 
visible. Monte Carlo numerical simulation of dust emission indi- 
cates that coma and tail were optically dominated by dust grains 
greater than 10 micron. Mass-loss rate has been estimated to be 
of the order of 0.2 kg s _1 . Dust emission probably began at least 
few months before the discovery of the activity (October 2005). 
This estimation of dust emission rate cannot be directly com- 
pared with the results of the 5 Cases shown before, because 
we have no idea if the activity is confined at the bottom of a 
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crater or if it is instead is widespread on large areas of the sur- 
face. In order to give some numbers, let us make the assumption 
that only the bottom of a recently excavated crater at the subso- 
lar point is the source of the observed activity (the walls would 
not contain volatile material). The surface area of P/Read is ap- 
proximately l.lxl0 6 m 2 . If we assume a dust production rate, 
at the equator, of 8.8xlCT 6 kg s~'irT 2 (figure[2]), in order to ob- 
tain the observed dust emission rate we need an active area of 
approximately 2.3.xl0 4 m 2 . The evaluation, in turn, gives a very 
rough estimation of the depth of the crater, if we assume a bowl- 
like shape. The implied crater diameter, for a planar surface, is 
~ 170 m. Assuming a depth-to-diameter ratio of 0.2 (Melosh 
1989), the corresponding crater depth is ~ 34 mThis number 
cannot be interpreted as a precise measure of the thickness of 
the original mantle of P/Read (could have been thinner at the 
impact point), but is only useful to give an idea of the orders of 
magnitude involved. This rough computation is performed for a 
crater close to the subsolar point: in case of an impact at a higher 
latitude, all the numbers should be scaled correspondingly. 
The crater size estimate given above is a lower limit. It would 
apply only for inert mantle thickness much lower that the esti- 
mated crater depth. However, we expect that MBCs have mantles 
thickness of the order of tens of meters because thinner mantles 
would have not have been able to prevent the complete sublima- 
tion of volatiles for MBC-sized objects. For instance, in the case 
that a thickness of half of excavated depth consist of inert man- 
tle, the same active area is obtained for a crater diameter about 
V(2) times larger. 

Looking to the plots in figures Q] [2] [5] and [6] the hypothesis of 
an impact leaving a mantled bottom appears much less compat- 
ible with observation, unless we consider a surviving dust layer 
much thinner that in the Cases 2, 3, 4 and 5: the activity level is 
so faint that it would be hardly detected from ground. 

3. How frequent are the impacts? 

In order to address if cratering events are a viable process to ex- 
plain the observed activity and number of MBCs, it is important 
to estimate not only the size of the crater required to sustain the 
observed activity, but also the frequency of its formation. 
The estimate of the projectile size needed to form a crater 
of 170 m has been performed considering two scaling laws, 
namely for porous material and cohesive soils (Holsapple and 
Housen 2007). The material strength has been set to 2xl0 4 and 
10 7 dyne/cm 2 . We assumed a target and projectile density of 
1000 kg/m 3 and 2600 kg/m 3 respectively. Moreover, the com- 
puted P/Read intrinsic average impact probability within the 
Main Belt is about 3.48xl0 18 km 2 yr~l and the average impact 
velocity is about 3.7 km s . These values are computed accord- 
ing to the Bottke et al. ([2005 ) Main Belt population, and follow- 
ing the procedure described in the work of Marchi et al. ((2010). 
The derived projectile sizes are ~ 8, 18.5 m for porous and co- 
hesive material, respectively. 

The computed time for having a collision with a projectile > 8 m 
on a body with the size of P/Read (0.6 km) is 5.1xl0 7 yr. The es- 
timated total number of objects larger than 0.6 km is 2.9xl0 6 , 
therefore such a collision is expected to occur on average every 
~ 17 yr everywhere in the Main Belt. This time become ~ 147 yr 
for projectiles > 18.5 m. 

Of course these numbers have to be taken as a rough indication, 
given the uncertainties in the Main Belt size distribution at such 
small size and in the target material parameters. Note that the 
timescales above also depend on the target size. For a body hav- 
ing the size of Elst-Pizarro (5 km), it is estimated that the same 



impacts occur every ~ 10 and ~ 100 yr, respectively for projec- 
tiles of 8 and 18.5 m. 



4. Discussion 

When running a simulation, it would be impossible to test all the 
possible combination of input parameters. What is usually done 
is, once found a meaningful combination, to keep fixed values 
for most of the parameters and vary only one or two at a time, 
in order to test the influence of the changes on the results. The 
variable parameters are the ones considered to be most influen- 
tial on the properties or phenomena that are the object of the 
simulation. In this paper, an initial composition for P/2005 Ul 
(Read) has been assumed, and we have varied, to build the 5 
Cases that have been studied, the parameters defining the thick- 
ness and composition of the dust mantle. As for the possible 
effects of changing the value attributed to the remaining initial 
parameters, it is here discussed on the basis of our experience. 
A lower value of dust/ice, (for example, a value of 1 , usually at- 
tributed to classical comets) would give as a result a lower dust 
flux, and consequently a slower erosion of the surface. Changing 
the thermal conductivity, porosity, density, pore radius would 
have the effect of slowing or quickening the evolution of the 
model body. The thickness of the surviving mantle in the Cases 
2, 3, 4 and 5 has been chosen in a quite arbitrary way. It is clear, 
anyway, that a larger thickness would tend to completely quench 
any possible activity. The mantle thickness necessary to com- 
pletely block the underneath sublimation strongly depends on 
the average porosity, thermal conductivity and the other physi- 
cal properties of the dust layer. 

The issue of the initial temperature attributed to the model body 
deserves some more discussion. Surface temperature is mainly 
determined by the solar input, and is spanning, along the or- 
bit followed by the model body, the range 147 - 190 K (subso- 
lar point, fast rotator approximation). Obviously, nobody knows 
what really is the average internal temperature (below diurnal 
and seasonal skin depth) of such a body, but the assumption has 
been made that it should be around 130 K, a value at which wa- 
ter ice can survive. When the insulating mantle is disrupted, the 
internal temperature is steadily rising till an equilibrium value is 
reached. A different value of the initial temperature would result 
in a slowing or quickening of this process. 
In this paper the possible effect of a single impact has been simu- 
lated. During the long time spent orbiting around the Sun, many 
small impactors hit the surface of bodies in the Main Belt. Even 
if these impacts are not able to disrupt the insulating mantle, it 
is probable that they are producing some changes on the surface. 
The effect of the impact itself, and of the released heat, could be 
that, for example, of locally compacting a porous surface mantle 
and locally changing its properties. 

The results obtained for P/ 2005 Ul (Read) can be easily applied 
to the other MBCs known until now, because they are all revolv- 
ing on similar orbits, with the only possible exception of P/2008 
Rl (Garradd), that is coming closer to the Sun at perihelion. For 
this object, solar input is slightly higher and the results obtained 
for P/2005 Ul (Read) should be consequently scaled. 
From the results of these simulation (in particular of Cases 2 to 
4) , it is also possible to understand how icy material buried deep 
enough under an insulating mantle could have survived for such 
a long span of time. A thorough analysis of this process will be 
the subject of a future paper. As for the origin of these ice-rich 
layers still surviving in bodies orbiting well inside the snow line 
border, it must probably be sought in the chaotic processes that 



M. T. Capria et al.: The activity of Main Belt comets 



9 



were shaping the solar system at the beginning of its history, 
bringing bodies away from their formation places. 

5. Summary and conclusions 

Assuming that the newly discovered Main Belt Comets are icy 
bodies on which a stable cometary activity has been induced by 
a recent impact that excavated a crater in a dust mantle, two dif- 
ferent scenarios have been simulated. In the first one, it has been 
assumed that ice-rich layers, remained buried under an insulat- 
ing dust layer for a very long period of time, are now directly 
exposed to Sun heating. In the second one, it has been assumed 
that the bottom of the newly excavated crater is still covered by 
a layer of dust much thinner than before the impact, so that the 
solar input is now able to reach ice-rich layers. The characteris- 
tics attributed to the model body are those of P/2005 Ul (Read). 
It has been demonstrated that, given the existence of buried ice, 
the activity observed in the MBCs can be explained as cometary 
activity, that means dust emission sustained by ice sublimation. 
In the first scenario, a stable cometary activity can be obtained, 
compatible with observations. Making a rough estimate, a crater 
at the subsolar point with a diameter of 170 meters would be 
compatible with the observed dust production rate. The activity 
would be stable, but the water ice is quickly consumed, depend- 
ing on the cometocentric latitude of the impact point. At the sub- 
solar point, the surface erosion could reach two meters per orbit. 
When, as assumed in the second scenario, a dust layer still re- 
mains, covering the bottom of the crater, a stable activity is still 
possible, providing that the heat of the Sun is able to reach ice- 
reach layers underneath. The activity is much fainter and long 
lasting than in the preceding case, depending on the latitude of 
the impact and the physical properties of the dust mantle. This 
kind of activity is probably not compatible with the activity ob- 
served on P/2005 Ul (Read). 

The computed time for having a collision with a projectile > 8 m 
on a body with the size of P/Read (0.6 km) is 5.1xl0 7 yr. The es- 
timated total number of objects larger than 0.6 km is 2.9xl0 6 , 
therefore such a collision is expected to occur on average ev- 
ery ~ 17 yr everywhere in the Main Belt. This time become 
~ 147 yr for projectiles > 18.5 m. For a body having the size of 
Elst-Pizarro (5 km), it is estimated that the same impacts occur 
every ~ 10 and ~ 100 yr, respectively for projectiles of 8 and 
18.5 m. 

There are no clues as how many other MBCs could exists in the 
Main Belt. Surely the great majority of bodies still keeping some 
ice in their interior are presently inactive, or their activity is so 
faint that it cannot be detected, at least from Earth-based obser- 
vation. 
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